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We study dynamical systems that switch between two different vector fields depend- 
ing on a discrete variable and with a delay. When the delay reaches a problem- 
dependent critical value so-called event collisions occur. This paper classifies and 
analyzes event collisions, a special type of discontinuity induced bifurcations, for 
periodic orbits. Our focus is on event collisions of symmetric periodic orbits in 
systems with full reflection symmetry, a symmetry that is prevalent in applications. 
We derive an implicit expression for the Poincare map near the colliding periodic 
orbit. The Poincare map is piecewise smooth, finite-dimensional, and changes the 
dimension of its image at the collision. In the second part of the paper we apply 
this general result to the class of unstable linear single-degree-of-freedom oscillators 
where we detect and continue numerically collisions of invariant tori. Moreover, we 
observe that attracting closed invariant polygons emerge at the torus collision. 
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1 Introduction 

Relay control systems, which can be regarded as hybrid systems, are applied in many different 
areas of engineering applications. Relay feedback control might involve control of stationary 
processes in industry as well as control of moving objects, for instance in flight control. Hence, 
it is not surprising that in recent years much research effort has been spent on investigations of 
the dynamics of relay systems characterized by idealized off/on relay control, see for instance 
[11, 16] among other works. However, relay systems often feature intrinsic hysteretic behavior 
[22] as well as delay in the control input. For instance, if we control a plant using simple on/off 
relay feedback control via a computer network, our control input may be delayed, due to the 
traffic along the network. 
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In the current paper we study the dynamics of relay control systems with hysteretic behavior 
and delay in the control input. We wish to understand if and how the interplay between the 
discontinuous nonlinearities (switching events), hysteresis and delay in the switching function 
can lead to new bifurcations. We expect so-called event-collisions to occur as observed in [7]. 
We show that they cause a change in the dimension of the phase space that describes the local 
dynamics of the system, which leads to interesting dynamical phenomena such as closed invariant 
polygons and collisions of smooth invariant tori. 

Let us explain the problem for the classical prototype example of an oscillator subject to a 
relay switch: 

x-\- tl,x-\-(0^x = u. (1) 

This example will be studied in detail in Section 5 and in Section 6. In (1) j = is the 

continuous variable and u provides discrete feedback in the form of a relay. That is, u is ruled by 
the switching law: 

M if±h'y{t-z)>e 
u{t) = < (2) 
\]im,/tu{s) if h^y{t-z) e{-e,e). 

Roughly speaking, the switching law (2) means that we set m to — 1 whenever we observe that 
h^y > £ and we set u to +1 whenever we observe that h^y < —e (negative feedback). The region 
< £} provides a hysteresis in the relay (a buffer between subsequent switchings from u 
to —u): we leave u at its value when we observe that < e. The lines {y : h^y = ±£} are 
called the switching lines or, more generally for higher dimensional y and nonlinear switching 
functions h, switching manifolds. An ordinary differential equation such as (1) with a switching 
law of the form (2) is one of the simplest forms of hybrid dynamical systems [21]. The twist in 
problem (1), (2) comes from the presence of a delay T in the observation (or the switching) in (2). 
Figure 1 illustrates how a trajectory of (1) might look like. The presence of a small delay T is 
often equivalent to a perturbation of the switching law (only if £ is positive, see, for example, [7]). 
For large delays T the dynamics shows an enormous amount of complexity [4, 5, 7, 15] because 
the number of switching times in the interval [t — r,t] is effectively included into the dimension 
of the system which then becomes larger than the dimension of the physical space [27] . 

This paper studys the transition from 'small' to 'large' delay. This transition is characterized by 
event collisions, which are discontinuity induced bifurcations [7]. The notion of a discontinuity 
induced bifurcation can be made precise for periodic orbits of (1), (2). The possible types of 
codimension one discontinuity induced bifurcations for generic systems with delayed switching 
have been classified and analyzed in [27] (without hysteresis): 

A. (grazing) the periodic orbit touches the switching manifold quadratically without intersect- 
ing it, 

B. (corner/event collision) the continuous variable y is on the switching manifold at two 
times with difference equal to the delay t: say, at time t and time t — r. The name is due to 
the fact that the periodic orbit typically has a 'comer' at t because it switches the vector 
field (u I— > —u) at t. 

Most of the event collisions reported in [7] violate one of the secondary genericity conditions in 
[27], namely that only one corner collides. The reason behind this unexpected degeneracy of the 
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Figure 1 : Sketch of qualitative behavior of hybrid systems with delay. Initial history segment 
of (x,i) is assumed to be constant, initial value of m is +1. The solid lines are the 
switching manifolds h = ±e. At the times tj the trajectory crosses the switching line, at 
the times tj + T it switches the control u to —u. 

phenomena observed in [7] is that the oscillator system (1), (2) (which is the subject of study in 
[7]) is symmetric with respect to reflection at the origin: 

y 1-^ —y, u ^ —u. (3) 

The periodic orbits of primary interest inherit this symmetry. If the control switches from u to —u 
at time t it switches from —u to u one half -period later (at t + p/2ii pi?, the period). Whenever 
the comer at t is on the switching manifold the opposite corner {&tt + p/2) is on the switching 
manifold as well due to reflection symmetry. This is in contrast to the secondary genericity 
conditions in [27]. In fact, all systems studied in [4, 5, 7, 15] have a reflection symmetry because 
they are piecewise linear, which renders the generic theory inapplicable to symmetric periodic 
orbits in these examples of practical interest. 

We close this gap in the general context, that is, for 3^ € M" and nonlinear right-hand-sides and 
switching laws. The first main result of the paper is an implicit expression for the local return 
map for symmetric periodic orbits close to a simultaneous collision of two corners. The simplest 
example for this scenario is a symetric periodic orbit of period p = 2z for which 3^(0) and y{r) lie 
both on the switching manifold. Due to the reflection symmetry this event has codimension one. 
The return map near a colliding orbit is a piecewise smooth map 7^ that has a (?i — 1) -dimensional 
image in one subdomain of its phase space and a n-dimensional image in the other subdomain. 
This type of map has been studied extensively in the context of grazing bifurcations in Filippov 
systems (systems where e = T = 0; see the textbook [9] and the recent review [18]). In this sense 
our paper provides a Hnk between the theory of systems with delayed switching and hysteresis to 
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the theory of low-dimensional piecewise smooth maps as treated in [9]. 

The second aspect of the paper is a study of the concrete class of oscillators (1), (2) near event 
collisions of symmetric periodic orbits. We pay special attention to parameter regimes where the 
event colhsion leads to quasi-periodicity and, thus, colhsions of invariant tori and the appearance 
of closed invariant polygons. None of the phenomena described in the sections 5 and 6 are present 
in [7] because the authors restricted their study to the special case of pure position feedback 
/i^ = (1,0) in the switching law (2). 

The first part of our paper adopts a 'local' approach, which is slightly different from the studies 
[4, 5, 7, 14, 15]. We do not aim to classify the dynamics of a particular class of systems as 
completely as possible. Instead, we develop a local bifurcation theory, considering a general 
system with ?i-dimensional physical space and assuming that it has a periodic orbit {y^, , ) . Then 
we study the dynamics near the periodic orbit deriving conditions for the presence of bifurcations 
of a certain codimension. In this way our results will be more general than studies of specific 
classes of systems but all statements are valid only locally. The consideration of only two vector 
fields and and a binary switch is not really a restriction when one studies the local dynamics 
near a particular periodic orbit. The second part of the paper then demonstrates how this local 
bifurcation theory can be useful in combination with numerical continuation to understand the 
dynamics of the class of oscillators (1), (2). 

The paper is organized as follows. In Section 2 we introduce the basic notation and collect 
some fundamental facts about the forward evolution defined by general systems with delayed 
switching and hysteresis. We also point out differences to the neighboring cases of zero hysteresis 
or zero delay. Section 3 specifies conditions on periodic orbits which guarantee that a local 
return map exists and that this map is finite-dimensional and, generically, smooth. Section 4 
classifies and analyzes the codimension one discontinuity induced bifurcations of periodic orbits 
for systems with reflection symmetry. Section 5 studys the single-degree-of-freedom oscillator 
(1), (2) near an event collision, classifying the local bifurcations of symmetric periodic orbits. 
In Section 6 we unfold the collision of a Neimark-Sacker bifurcation with the switching line, 
a codimension two event involving collisions of closed invariant curves and closed invariant 
polygons. The appendix contains the proofs of all lemmas and the technical details of some 
constructions. 



2 Fundamental properties of the evolution 

We consider general hybrid dynamical systems with delay of the following form: 



m 

u{t) 



m)At)) 

-1 if h{y{t - t)) > £, or 

if h{y{t - t)) G (-£,£) and u^{t) = -1, 
1 if h{y{t - t)) < -£, or 

if h{y{t - t)) G (-£,£) and U-{t) = 1 



(4) 



(5) 



where £ is positive. The continuous variable y is ?i-dimensional and the discrete variable u is, 
for simplicity of presentation, binary, controlling the switching between the two vector fields 7j_ 
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given by y = /(j, ±1). In the definition of u{t) in (5), u^{t) is defined as 

:= limu{s). (6) 

We assume continuous differentiability for the functions /(•, ±1) : M" M" and /j : M" M in the 
right-hand-side of (4), (5) with respect to the argument y (and, possibly, parameters). Furthermore, 
we assume that the gradient h' is non-zero everywhere and that / and h have uniform Lipschitz 
constants (using a prime for the derivative with respect to y): 

11/11 < Lm^x, and \\h'\\ < //max- 

Due to the delay in the argument of y in the switching decision (5) the phase space of (4), (5) is 
infinite-dimensional. An appropriate initial value fory is the history segment j([—T,0]) [12, 29]. 
Thus, the phase space of (4), (5) is C([— 0,O];M") x { — 1,1} where is an upper bound for 
the delay z (we will vary T as a bifurcation parameter in the sections 5 and 6). The notation 
C([— 0,O];M") refers to the space of continuous functions on the interval [—0,0] with values in 
W. 

We clarify in Appendix A in which sense system (4), (5) constitutes a dynamical system. 
In particular, we give a precise definition of the forward evolution E*{^,uq) from an arbi- 
trary initial condition {^q,uo) G C([— 0,O];M") x { — 1, 1} using the variation of constants for- 
mulation of (4), (5). The evolution E'{^o,uo) has a continuous infinite-dimensional compo- 
nent E'^{^Q,uo) G C([— 0,O];M") and a discrete component E'^{^o,uo) G { — 1,1}. The infinite- 
dimensional components are related to each other by a simple time shift: 

E'Mo,uom=E'+\MiO) 

with a common headpoint trajectory EI.{^o,uq){0) G C([0,f£];M") for all 6 G [—0,0] and t G 
[O,?^]; see [12]. Thus, ^^(i^OiMo) is continuous in f for all f > 0. However, in general £^,(1^01 mo) 
(for a fixed t) does not depend continuously on the component of its initial value. 

Lemma 1 summarizes two basic facts about the evolution E'{^,u) that are proved in Ap- 
pendix A. 

Lemma 1 (Fundamental properties of evolution) 

Let E' {t,Q,uo) be a trajectory of the dynamical system (4), (5) on a bounded interval [tQ,tE]- Then 
the following holds. 

1. The discrete component E'j{^o,uo) changes its value (switches) only finitely many times in 
[tcijE]- Thus, Ec follows either or Y^for all but finitely many times t G [foi^fi]- 

2. Ift() > T then the number n^ of switchings ofE^ is bounded by 

ns <! + [?£— to] • Lmax-^^maxJmax/ (2£) 

where y^ax the maximum of [[^^(i^cMo)!,^^,,,,,^] || and L^ax and //max are the Lipschitz 
constants of f and h, respectively. 
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The statements in Lemma 1 are subtly dependent on the presence of hysteresis and delay, which 
can be seen from the fact that statement 1 is not true in general if e = 0. See [27] for details about 
systems with switches that have delay but no hysteresis and [4, 5] for studies of a piecewise linear 
oscillator with rather intricate behavior. 

Lemma 1 implies that the headpoint trajectory El.{^o,uo){0) is differentiable with respect to 
time t in any bounded interval [?o,?£], following one of the flows Y±, except in a finite number of 
times. 

3 Periodic orbits 

If the relay is used to control a linearly unstable system, such as the oscillator in (1) with negative 
damping i^, we cannot expect to have stable equilibria. In this case the simplest possible long- 
time behavior of (1) is periodic motion. We focus on the dynamics near periodic orbits for the 
remainder of the paper. This section introduces the necessary notation and collects some basic 
facts about the behavior of the general relay system (4), (5) near a periodic orbit. 

Definition 2 (Crossing time) Let L = (yi,{t),u^{t)) be a periodic solution of the general relay 
system (4), (5), that is, y.f{t + p) = y*{t), u-^^it + p) = u^{t)for allt >0 and some period p > 0. 
We call a time t crossing time ofL if\h{yif{t)) \ =£ and M*^_(f + t) = lim^^f m*(5+ t) ^u^{t + 

Equivalently, we could say that ? is a crossing time if h{y^{t)) = £M*._(f + t) = £\mis/t u^{s + x). 
Time T after a crossing time will switch. A periodic orbit L can have only finitely many 
crossing times tk{k= 1, . . . ,m) per period. Thus, L is differentiable, following one of the flows 
Y±, in all times except, possibly, mtk + T + jp {k = 1 , . . . , ni, 7 G Z). We assume that the period 
p is larger than the delay T (without loss of generality because p does not have to be the minimal 
period). 

First, we establish when the evolution is continuous with respect to its initial value in a periodic 
orbit. We call the condition for continuity weak transversality. As the name suggests it excludes 
that the periodic orbit touches the switching manifold without crossing it (but does not require 
positive speed of crossing, thus, weak transversality). If this condition is satisfied then it makes 
sense to define a local return (or Poincare) map along the orbit 

Definition 3 (Weak transversality) Let L = {y^{t),u^{t)) be a periodic solution of period p of 
the general relay system (4), (5). We say that L satisfies weak transversality if\h(y^,(t)) \ is locally 
strictly monotone increasing near all crossing times t\,...,tm^ ^^p] ofL. 

Weak transversality permits that lim.,/-; 3)* {s) / lim^^, {s) and each of the limits may be tan- 
gential to the switching manifold {h = ibe} in any crossing time t of L. However, it enforces that 
the periodic orbit L cannot just touch the switching manifold quadratically, and that the number 
m of crossing times is even. Weak ti^ansversality is generically satisfied for a periodic orbit. 

Lemma 4 (Continuity) If a periodic orbit L = (3^* is weakly transversal then the 
continuous component El(^,u) of the evolution E is continuous with respect to £, in {£,,u) = 
{y.^{-),u^{0))forallt > 0. 
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(See Appendix A for proof.) Due to this continuity it makes sense to define a local return map 
(also called Poincare map) along a weakly transversal periodic orbit L = {y^{t),u^{t)) to a local 
cross section y. Assume (without loss of generality) that neither nor — T is a crossing time 
of L and that follows F+ near time (thus, m*(0) = 1). Then is differentiable in 0. Let 
L = {y€W: y*{Of[y-y*{0)] = 0} be the hyperplane in M" orthogonal to %{0) in 3;*(0). We 
choose as Poincare section 

^ = {(§,M)GC([-0,O];/?")x{-l,l}:^(O)er, u = \}. (7) 

It is not actually necessary to choose £ orthogonal to j*(0). Any cross-section £ which is 
transversal to 3^* at f = is admissible. Let — T < 1 < . . . < t^ ^ < < < . . . < f*,,,, < p — T 
be the crossing times of L in [—T,p — r]. The following lemma states that the local return map 
to is in fact a map in £ x M'^ (a space of dimension n — l+ jj.). The notation U (v) refers to a 
(sufficiently small) neighborhood of a vector or number v in Lemma 5 and throughout the paper. 

Lemma 5 (Finite-dimensional Poincare map) 

There exist a 5 > and neighborhoods Ui{y{-)) C Ujiyi)) C C([-0,O];/?") and U{p) C M 
such that the following holds: All initial conditions {E, ,u) ^ = ^ r\ \U\{y{-)) x {1}] have 
a unique return time T{t,) £ U {p) to 5^2 = y P\ \U2{y{-)) x {I}]- For any initial condition in 
(i§,m) G Ui{y{-)) X {1} there exist unique times t[ < . . . < t^ in (— T,0) such that 

tj = mm{t £[t,j-5,t,j + 5]:\h{^{t))\=e}, (j=l,...,^). (8) 

The local return map 

depends only o?i (1^ (0) , . . . , f^) G £ x 

We have omitted the discrete component of the return map from P because it is always + 1 . 
An equivalent definition of tj is tj = min{f G — 5,t^ j + 5] : h{^{t)) = euj} where uj = 
u*-{t^j + t). The precise dependence of P on {^{0),ti, ...,t^) is given in Appendix A. 

Corollary 6 (Smooth Poincare map for generic periodic orbits) Assume that all crossing times 
t*,k (k = ■ ■ ,tn) of the periodic orbit L= (y^,u^) satisfy the following two conditions: 

1. (no collision) the time t^ j- — z is a not a crossing time ofL, and 

2. (smooth transversality) /z'(3'*(?*,/t))j* / 0. 

Then the Poincare map P depends smoothly on the coordinates (i^(0),fi, . . . ,f^) G £ x M.^ of 
,u) £ y (ti, . . . ,t^ as defined by (8)j. 

The two conditions of Corollary 6 guarantee that whenever {t) changes its value then y^ {t — t) 
follows one of the two flows Y± in a neighborhood of f — T and {t — t) points transversally 
through the switching manifold {h = ib£}. Generically, these two conditions are satisfied, which 
implies that, generically, the Poincai^e map of a periodic orbit is smooth. Condition 2 is more 
restrictive than weak transversality as it requires differentiability of and a non-zero time 
derivative of /i(j*(-)) in 
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Corollary 6 implies that the dynamics and possible bifurcations near aperiodic orbit (y*,^*) 
of the general delayed relay system are described by the theory for low-dimensional smooth 
maps whenever the conditions 1 and 2 are satisfied. See, for example, [20] for a comprehensive 
textbook on bifurcation theory for smooth systems. 

Definition 7 (Slowly oscillating periodic orbit) 

We call a periodic orbit L = (y^,u^) slowly oscillating if the distance between subsequent crossing 
times tp- ofL is always greater than the delay T. 

If we choose the Poincare section y appropriately then the local return map of a slowly oscillating 
periodic orbit L satisfying the genericity conditions of Corollary 6 is simply a return map to Z, the 
hyperplane of W in the definition of J^. This has been observed in [7] for the oscillator (1), (2). 

4 Discontinuity induced bifurcations 

Let us assume that the right-hand- side / and the switching function h depend on an additional 
parameter A. What happens to the dynamics near a periodic orbit under variation of A (or, 
alternatively, the delay t)? The previous section has established that, as long as the conditions 
of Corollary 6 are satisfied, we should expect standard bifurcation scenarios such as period 
doubling, saddle-node or Neimark-S acker bifurcations (see [20] for a classification). However, 
when varying the parameter A we can also achieve that any of the conditions 1 and 2 fails at 
special parameter values. We call these events discontinuity induced bifurcations. 

For compactness of presentation we assume that for A < Ao the periodic orbit L = (3^* , ) is 
slowly oscillating, that is, the distance between subsequent crossing times 4. is always greater 
than the delay T for A < Aq. 

4.1 Generic bifurcations 

Generic grazing If condition 2 of Corollary 6 is violated at A = Aq the orbit L grazes the 
switching manifold tangentially at a crossing time . Let us assume that y^ follows in t^ 
without loss of generality. That is, 

0= ^/i(y*(0,^)|r=r, = /z' ),A()) =h'{y^{t^),?H))f{y^{t^),+\,A^). 
Generically, one can expect that 

/ jj^h{y^{t),A^)\,=t, = /z' Ao)j;(f*) +/j"(y,(f,),Ao) [Mt*)f' (9) 

which means that the periodic orbit L touches the switching manifold quadratically and not to a 
higher order. However, under condition (9) L is not weakly transversal. In general, we cannot 
expect that the evolution is continuous in L. Thus, trajectories arbitrarily close to L leave a fixed 
neighborhood of 3^* in a finite time (typically one period). Consequently, generic grazing of the 
periodic orbit L cannot be described by the approach of local bifurcation theory adopted in this 
paper. 
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Generic corner collision If condition 2 of Corollary 6 is violated at A = Ao the orbit L switches 
between the two vector fields exactly at a crossing time fi. That is, 

/ limM*(s) = -{h)- 

Let us denote f\ = f{y^{t\)^u.^^^{t\),%Q) and /2 = /(j*(fi),M*(fi),Ao)- Generically, we can 
expect that 

CI. all other crossing times t2,---,tm of L do not collide, that is, u^{tic) = u^-{tk) for k = 
2,.. . ,m, and 

C2. the left-sided tangent /i and the right-sided tangent fj to y^{-) in are transversal to the 
switching manifold. More precisely, 

q:=h'iy,{h),Xo)fi- h'{y, (fi ) , Ao) /z / 0. (10) 



If ^ < in condition (10) the periodic orbit L has a comer at ti and this comer touches the 
switching manifold from one side at the crossing time ti. Thus, L is not weakly transversal to 
the switching manifold in its crossing time ti. In this case we cannot expect that the evolution 
is continuous in L. Consequently, comer collisions with ^ < cannot be described using local 
bifurcation theory, either. An explicit expression for the retum map (which is a piecewise 
asymptotically linear {n — 1) -dimensional map) under the assumptions CI and C2 for the case 
^ > has been derived in [27] (as case (a) in Appendix D of [27]). 



4.2 Corner collision with reflection symmetry 

Often the practically relevant examples have special symmetries which enforce that if one corner 
of a symmetric periodic orbit collides at A = Ao then other comers of the symmetric periodic orbit 
collide at Aq simultaneously. For example, the systems studied in [4, 5, 7, 15] and our prototype 
oscillator (1) are all piecewise affine: 

f(y, u) =Ay + bu 

r (11) 

h{y)=h'y. 

Thus, they have a full reflection (Z2) symmetry 

/(y,M,A) = -/(-J,-M,A) 

h{y,X) = -h{-y,X). 

This means that, even though condition C 1 should be generically satisfied for a colliding periodic 
orbit, restricting to the generic case disregards many practically relevant systems. 

The Z2 symmetry (12) typically gives rise to a symmetric periodic orbit L = {y^,u^) satisfying 
3'*(f — T) = —y>f{t) and u^{t — T) = —u^{t) for the half -period T and all times t. A comer 
collision of L for a crossing time f at a special parameter A = Ao automatically induces a corner 
collision for the crossing time t — T, a. scenario that is not covered by the generic bifurcation 
scenarios listed in Section 4. 1 . 
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Let us assume that system (4), (5) has full reflection symmetry (12) and a symmetric periodic 
orbit L = (j* , ) of half-period T that experiences a corner collision at the parameter A = Aq for 
crossing time and, enforced by symmetry, for crossing time T . For compactness of presentation 
let us assume that and T are the only crossing times of L. This implies that the delay T equals 
the half -period T and that switches between +1 and —1 at the crossing times and T . Without 
loss of generality L consists of the two segments 

y,{[0,T]) =f|''^13;,(0), M*([0,r)) =+1 

y,{[T,2T]) = yi^'^V^Cr) = -i-f ^^j*(0), M.([r,2r)) = -i 

Moreover, h{y^{0),XQ) = e and /j(j*(r),Ao) = — £■ Note that the colliding orbit L is always 
following the 'wrong' flow. That is, the orbit is identical in shape to the periodic solution with 
positive feedback (+1 and —1 interchanged in (5)) and zero delay. In addition we assume that 
condition C2 is satisfied for the collision time f = (and, by symmetry, for t = T) with ^ > 0. 
We call this condition strict transversality because it is stronger than the weak transversality 
introduced in Section 3: 

q:=Kfi-Kf2>0 (13) 

where = h'{y, (0) , h) , fx = f{y* (0) , - 1 , Aq) and /2 = f{y, (0) , 1 , A(,). 

This strict transversality guarantees that the evolution of the continuous component ^'((^ , u) is 
continuous for ^ =y*{s + -) for all G M. We choose a cross section £ for the Poincare map at 
f=A(0<A^l) and orthogonal to /2: L = {y : f2 [y — y*{A)] =0}. If A is sufficiently small 
then £ is transversal to 3^^, in f = A because y^ (A) = /i + 0(A) . 

Using the cross-section £ in the definition of the Poincare map P, Lemma 5 states that the 
return map P in its domain of definition 

^1 = [/(y, (A + •)) n G C( [-0, 0] ; M") : ^ (0) e £} 

depends for an initial value ^ € J/^i only on the headpoint ^ (0) G £ and the time ?i ) = min{f G 
[— 2A,0] : h{^{t)) = e}. The time ti{^) is the time when ^ crosses the switching manifold 
{h = e} for the first time in [— 2A,0]. This time ti{^) exists if U{y^{A + •)) is sufficiently small. 
Effectively, the map P depends only on n coordinates. 

The following lemma simplifies the representation of the Poincare map P to a map from 
[/(3;*(0))backto[/(3;*(0)). 

Lemma 8 (Return map near collision) All elements ^ of the image rgP of the return map P 
have the form 



y^eCv)]^ /f^G [-0,-0(3;)], 



^w=rw viv::"?:. (14) 



where y & U (3'*(0)) and d{y) is implicitly defined by the condition Y^^^\ G £. The coordinates 
{t, (0),fi(i§)), which are necessary for the definition ofP, are uniquely defined by y. 
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(See Appendix B for the proof.) We have omitted the discrete variable u as an argument of P 
because it is always +1 at the cross-section. The point y = ^{—d{y)) is the point where the 
continuous component G rgP switches from following F_ to following F+. The time d{y) is 
the time that has elapsed between the switching and the intersection of the headpoint with Z. For 
a sufficiently small ?7(y*(A + •) we can assume that d{y) is in [0,2A]. Lemma 8 states that the 
return map P, restricted to its image TgP, can be described as a map from the switching point of 
the initial value to the switching point of its image under P. 

Exploiting the representation (14) we can describe the dynamics of the Poincare map P by a 
map m mapping from U {y^ (0) ) C W back to U (y* (0) ) . 

Theorem 9 (Reduced return map near collision) Let X he sufficiently close to Aq. The return 
map m for elements of the image rgP of the Poincare map P is given by m = F oF where 
F : [/(j*(0)) ^ U{y40)) is defined by 

F{y) = -rr^-^ (15) 
and t(y) € (—A, A) is the unique time such that 

e=h (y!_^^\) ifh{y) > £, 



■h(Yi^'\) ifh{y)<e. 



The expression (16) of the traveling time t{y) implies that F is continuous in U (3'*(0)) and smooth 
in each of its two subdomains 

D^ = U{yM)n{y:h{y)>e} and 

D+ = uiy,mn{y-Ky)<£} 

but, in general, its derivative has a discontinuity along the boundary Do between D_ and D+. 
The regularity of the implicit expression (16) for t{y) for all y G ?7 (3'*(0)) follows from the strict 
transversality (13) of the colliding periodic orbit L. Let us denote the two smooth parts of the 
map F by F^. and F : 

F+{y) = -F^+'+Wj where h (f+^^'^V) = £> and 

(IS) 

F_{y) = -F^+'-Wj where h (f^T^'V) = £■ 

That is, F\o^ = F^ and F\d_ = F . Due to the regularity of the definition of t±{y) in U (3'*(0)) 
the maps F± can both be extended to the whole domain D. The map F+ projects D nonlinearly 
onto the local submanifold (the delayed switching manifold) 

rgF+ = -Yl[{h = e}nU{yM)] = {y ^ U^yM) ■ KY+'[-y]) = £}, 

which has co-dimension one. The linearizations of F± with respect to y and the delay z in 

y=y^{0) are 



F±(3;,(O) + T];t + 0)=3;,(O)-A(t) 

g- = h'ofi = -h'^'fi 



g± 



ri + Bf2 



(19) 
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where = dyYl in 3;,(0), and /i = f{y, (0) , - 1 , Ao), /z = f{y, (0) , + 1 , Ao), and = /j'(3;, (0)) 
(as introduced before). The linearization of F+ projects rj by the linear projection / — /2/Jo/ {h'^fi) 
before propagating it with A (t), mirroring the dimension deficit of the image of the nonlinear 
map F+. 

The map F is a continuous piecewise smooth map in M" with a rank deficit in one half of the 
phase space. F is implicitly defined and nonlinear, even if the original relay system is piecewise 
linear of the form (11). Thus, Theorem 9 reduces the study of the dynamics near a colliding 
symmetric periodic orbit to the study of a low-dimensional piecewise smooth map in a similar 
fashion as for the generic case [27]. General bifurcation theory has been developed for piecewise 
affine maps in M", which carries over partially to the nonlinear case [2, 9, 10, 18, 23, 28, 30]. 



5 Single-degree-of-freedom oscillators 



In this section we demonstrate the use of the reduced maps derived in Theorem 9. We study the 
linear single-degree-of-freedom oscillator (1) subject to a delayed linear switch with hysteresis 
(2). We rescale time t and the variable x in (1) such that the equilibria of the flows Y± are at 
±1 and such that each flow rotates with frequency 1. Furthermore, we introduce the parameter 
a which is the tilting angle of the switching decision function h. This reduces system (1) to a 
system with four parameters, the damping ^, the delay T G (0,oo), the width 2e G (0,oo) of the 
relay hysteresis region, and the angle a e [— ;r/2, 7r/2] (negative feedback) of the normal vector 
to the switching lines: 

'2^--(l + ?')« 



A; + 2Ci + (l + C-)x 



u{t) 



-1 
1 

lim 



s/t 



ms] 



if x{t 
if x{t 
if x{t 



t) cosa - 
t) cosa - 
t) cosa - 



sin a > £ 
sin a < — £ 
sina G 



(20) 



We consider the case of an unstable focus (spiraling source) corresponding to < 0. In our 
numerical investigations we fix the damping to —0.1, which is a moderately expanding unstable 
spiral. We study the dynamics near the colliding periodic orbit and how it depends on the 
parameters of the control switch u (delay T, angle a and hysteresis width £). The affine flows Y± 
are given by 

Y!0=A{t)y±v{t) 

where 



A(t) = e cosf 
v{t) 



e ^'^sinf 



(^sinf 



1 + C' 

COS? + 



, and 



(1 + C')sinf 



and have equilibria at (±1,0)^. The condition for the existence of a colliding symmetric periodic 
orbit is 



3;icosa+j2sina = e for t > tt 
ji cosa + j2sina = — £ for t < tt 



(21) 
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Figure 2: Panel (a): visualization of a colliding periodic orbit of (1), (2). Panel (b): collision 
surface in the three-dimensional parameter space (T,a,£). Damping C, is fixed at —0.1. 
At the parameters on the surface there exists a symmetric periodic orbit L experiencing 
a comer collision. Beneath the surface the symmetric periodic orbit is a fixed point of 
the map in the subdomain D+, above the surface it is a fixed point of the map F in 
the subdomain D_ . 

where y = — [/ + A(t)]^'v(t) (thus, y = —Yj^y). For < 0, condition (21) ensures that the 
switching point y lies on the switching manifold and that /2 = /(j, 1 ) points out of the hysteresis 
region. See Figure 2(a) for a visualization of a symmetric periodic orbit satisfying the collision 
condition (21). Figure 2(b) shows the surface of parameters (t, a, e) where a colliding symmetric 
orbit L exists in system (20). Whenever one varies the system parameters along a path intersecting 
the surface transversally L undergoes a corner collision. Beneath the surface L is a fixed point of 
in D+, above the surface L is a fixed point of F in D . The change of stability of the fixed 
point of F at the corner collision is determined by the linearizations of F and F+. Figure 3 shows 
a map of the different eigenvalue configurations that can occur on the collision parameter surface 
for T G {7i,2n). The map F+ has a one-dimensional image. Thus, the linearization of F+ in its 
fixed point has only one non-zero eigenvalue. The horizontally hatched region in Figure 3 shows 
where this eigenvalue has modulus less than one. Within this region the fixed point of F+ is stable 
at collision. Consequently, the symmetric periodic orbit of (20) is linearly stable for parameters 
beneath the horizontally hatched region in Figure 3 of the collision surface in Figure 2(b). The 
stable region of F+ is bounded by a flip bifurcation (eigenvalue equals — 1, thin dotted in Figure 3) 
and a fold bifurcation (or saddle-node, eigenvalue equals 1, thin dashed in Figure 3). Parameter 
values on a bifurcation curve in Figure 3 correspond to codimension two events of the oscillator 
(20) because they also lie on the colhsion surface. This means that the symmetric periodic orbit 
of (20), as a fixed point of F, has a linearization with neutral stability in one of the subdomains 
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Figure 3: Bifurcations of the maps F+ and F projected onto the corner coUision surface ('look 
from above' onto the surface in Figure 2(b)). Damping is fixed at —0.1. Only the part 
of the collision surface with positive e is shown. Thick lines are standard bifurcations 
of F . Thin lines are standard bifurcations of F+. The cross is a coincidence of a 
Neimark-Sacker bifurcation of F with the comer collision at £ = 0.1. Its numerical 
unfolding is presented in Figure 4 and Figure 6. 

D± and, simultaneously, it is located on the boundary between D and D+. 

The map F has a two-dimensional image. Thus, the linearization in its fixed point has two 
potentially non-zero eigenvalues. The vertically hatched region in Figure 3 shows parameter 
values where both eigenvalues are inside the unit circle. In this region above the collision surface 
the oscillator has a symmetric periodic orbit that is linearly stable as a fixed point of F . The 
standard bifurcations of the fixed point of F are shown as thick lines (flip dotted, fold dashed, 
Neimark-Sacker hollow). Along the Neimark-Sacker bifurcation (also called torus bifurcation) 
curve a complex conjugate pair of eigenvalues is on the unit circle. All of the bifurcation curves 
of F correspond to codimension two events for the symmetric periodic orbit of the oscillator 
(20) because they occur simultaneously with the collision, lying on the collision surface in the 
three-dimensional parameter space shown in Figure 2(b). 

Remarks The flip bifurcation of F± corresponds to a symmetry breaking bifurcation of the 
original symmetric periodic orbit of the oscillator (20) because the return map along the full 
periodic orbit is the second iterate of F . 
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The line a = in the figures 2(b) and 3 corresponds to the case of pure position feedback 
studied in [7]. The linearizations of and F coincide for a = 0. More precisely, the expressions 
for g+ and g_ in (19) are identical. It has been observed in [7] that non-smooth phenomena 
cannot occur for symmetric periodic orbits at the event collision. 

The points (a, t) = (0, n) and (a, t) = (0,27r) in Figure 3 are highly degenerate. The periodic 
orbit does not intersect the switching line transversally at these parameter values, violating the 
strict transversality condition (13). The collision surface is singular for x = n. 

Figure 3 shows codimension two degeneracies of the linearization of F_ such as a concurrence 
of eigenvalues — 1 and +1 at PD-SN and strong resonances along the Neimark-Sacker bifurcation 
(double eigenvalue —1 at 1:2, eigenvalues exp(ib27r//3) at 1:3, eigenvalues it/ at 1:4, see [20] 
for an analysis and description). These points correspond to codimension three bifurcations of 
the symmetric periodic orbit. Similarly, all crossings of bifurcations of F and bifurcations of 
in Figure 3 correspond to codimension three bifurcations of the symmetric periodic orbit. Other 
bifurcations of higher codimension that can occur along the bifurcation curves shown in Figure 3 
involve the degeneracy of higher order terms in the normal form. These special points have been 
omitted from Figure 3. 

The interaction PD-SN between the flip and the fold of the fixed point of the map F_ appears 
degenerate in Figure 3 in the following sense. The flip curve crosses the fold curve transversally 
instead of touching it quadratically as one would expect in a generic parameter unfolding [19]. 
This is due to the projection of the bifurcation curves onto the collision surface. The set of all 
parameters in the (t, a, £)-space where the fixed point of F is a saddle-node forms a surface. 
This surface does not intersect the collision surface shown in Figure 2(b) transversally but touches 
it tangentially in the (dashed thick) fold bifurcation curves of Figure 3. The same applies to the 
fold of the fixed point of F^. 

6 Unfolding of Neimark-Sacker bifurcation collision 

The symmetric periodic orbit L of the oscillator (20) is stable near its comer collision in the 
hatched regions in Figure 3. Where the two hatched regions overlap L is stable for parameters on 
both sides of the collision surface in Figure 2(b). Of primary interest are the parameter regions 
near bifurcation curves bounding the region of stability of L. In these regions we can expect that 
other (possibly stable) invariant objects exist near L, which in turn collide with the boundary Do 
between the two subdomains D and D+ of the phase space. In order to understand the dynamics 
near one of the codimension two events one has to unfold it using two parameters. 

A systematic classification of possible unfoldings of codimension two bifurcations of piecewise 
smooth maps is not available in contrast to the situation for smooth systems [20] . Due to the 
impossibility of a general center manifold reduction it is difficult to derive general results. The 
references [2, 3, 9, 18] give a long list of possible cases but discuss unfoldings only for very few 
concrete examples. 

In this section we describe in detail the dynamics near the symmetric periodic orbit L near a 
collision of a Neimark-Sacker bifurcation of F with the boundary Dq. This case is only possible 
due to the increase of the dimension of the phase space from one (in D+) to two (in D ). In this 
sense it is the most characteristic feature of the symmetric corner collision of the oscillator (20). 
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Figure 4: Unfolding of colliding Neimark-S acker bifurcation in (t, a) -plane. Values of other 
parameters: £ = 0.1, = —0.1. All curves have been obtained by direct numerical 
continuation. The sketches on the right show qualitatively the dynamics in the different 
regions (a)-(d) (the dashed line is the switching line, attractors are solid, repellors 
hollow). Figure 5 shows the attracting invariant polygon at the parameter value marked 
by a cross. Figure 6 shows a parameter sweep of simulations along the gray arrow. 

It is also the most complex case due to the involvement of invariant curves. 

We fix £ = 0.1. The plane {(T,a,£) : £ = 0.1} intersects the collision surface of Figure 2(b) 
along a curve (not shown in Figure 3) and it intersects the Neimark-Sacker bifurcation on the 
collision surface (thick hollow curve) in a point NSC. This point NSC is marked by a cross in 
Figure 3 and corresponds to a set of parameters where the periodic orbit L is on the boundary 
Do, the linearization of F has a complex conjugate pair of eigenvalues on the unit circle, and 
the linearization of F+ is stable. We unfold this codimension two event using the parameter a 
(tilting angle of the switching line) and the delay T. Figure 4 shows the bifurcation diagram in 
the (t, a)-plane. The point NSC corresponds to the cross mark in Figure 3. The dynamics in the 
different regions is sketched in the insets next to the diagram in Figure 4. 

Let us explain the dynamics in the different regions near NSC (as sketched in the insets to the 
right of the diagram in Figure 4). In parameter region (a) L is in region D . L is a fixed point 
of F and it is linearly stable, its linearization having a pair of complex conjugate eigenvalues 
inside the unit circle. On the Neimark-Sacker bifurcation curve (hollow) this pair of eigenvalues 
is exactly on the unit circle. Thus, L changes its stability here. Moreover, the Neimark-Sacker 
bifurcation of L does not have a strong resonance (see Figure 3 where the strong resonances are 
marked). 

Consequently, a smooth closed invariant curve (invariant under F_) emerges from L at the 
Neimark-Sacker bifurcation near NSC. In this case the Neimark-Sacker bifurcation is supercritical, 
which implies that the emerging closed invariant curve is stable and exists in the region of linear 
instability of L (region (b) in Figure 4). As the diameter of the invariant curve grows it collides 
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Figure 5: Attracting piecewise smooth invariant curve. Parameter values: a = —0.44, T = 4.25, 
e = 0.1, = —0.1 (cross in Figure 4). Panel (a) shows the switching points of the 
trajectory of the oscillator (20) in the (jc,i-)-plane. The switching lines are drawn as 
solid lines. Panel (b) zooms into the neighborhood of the attractor of F. Panel (c) shows 
the circle map on the closed polygon. 

with the boundary Dq. This collision occurs along the solid curve in Figure 4. 

At this collision the smooth closed invariant curve disappears. In region (c) trajectories from 
all initial conditions close to L (except L itself) will eventually visit the region D+. Thus, they 
will eventually follow the map at least once getting projected onto rgF+, a one-dimensional 
manifold. Since the fixed point of is stable and in region D (thus, it is not a fixed point of 
F) all trajectories starting from points in D+ will eventually visit D . Hence, in region (c) all 
trajectories map back and forth between the two regions D and D+. 

This gives rise to invariant sets composed of finitely many smooth arcs that are images of a 
section of rgF+ under F_. Since F_ is approximately a rotation these arcs are rotations of rgT^^.. 
Their composition forms a piecewise smooth closed invariant curve consisting of finitely many 
arcs (a closed invariant polygon), which is located partially in D and partially in D+. This type 
of invariant set has been found and discussed for a piecewise linear map already in [30, 31], 
and for flows in [32, 33]. Figure 5 shows an example of such a polygon at the parameter values 
marked by a cross in Figure 4. Figure 5(a) shows the switching points of the trajectory of the 
original oscillator (20) (where the discrete component u changes its value). Figure 5(b) is a zoom 
into the vicinity of the upper switching region, showing the attractor of F . Figure 5(c) shows the 
map restricted to the invariant polygon parametrized by the angle (p of the point on the polygon 
relative to the the average of the attractor inside the curve. The map is clearly invertible giving 
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Figure 6: One-parameter pseudo bifurcation diagram for varying a. Values of other parameters: 
T = 4.2, £ = 0.1, = —0.1 (along the gray arrow of Figure 4). Envelope contains 
maximum and minimum of the iterates of F between 40 and 400 iterations. SPC, ICC 
and NS point out the intersection with the bifurcation curves in Figure 4. 

rise to either quasi-periodicity or a pair of periodic orbits (locking) on the polygon. According 
to [31] non-invertible maps on polygons (and, hence, in principle chaotic dynamics) are also 
possible but we did not encounter this phenomenon close to NSC. See also [30, 31] for a study of 
resonance and locking phenomena on the invariant polygons. 

If the parameters approach the collision curve of L (dashed curve in Figure 4) the invariant 
polygons shrink as the fixed points of F+ and F approach each other (and the boundary Do). 
Thus, in region (d) only the stable fixed point of F in D exists. 

Figure 6 shows the result of a parameter sweep along the gray arrow in Figure 4 showing 
the maximum and the minimum of the attractor (maxF^(3'o) and minF^(yo) for j = 40, . . . ,400, 
restarting such that 3^0 is the last iterate from the previous parameter). The figure gives evidence 
that the diameter of the invariant polygons grows linearly in the parameter a, starting from the 
point SPC where the corner collision of L occurs. During the sweep the polygon also changes 
its shape (for example, the number of comers or 'overshoot' of its edges). A detailed analysis 
how invariant polygons change under variation of the rotation and the linear expansion of the 
two-dimensional map F is carried out in [32]. The polygon attracts in finite time between the 
points SPC and ICC in Figure 6 such that the numerical result shown in the plot is very accurate 
for this region. The detection of the collision of the smooth invariant curve near ICC, the family of 
smooth invariant curves and the Neimark-Sacker bifurcation NS are, however, only inaccurately 
accessible by pure simulations due to the weak attraction in region Z)_ . 
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Figure 7: Family of closed invariant curves along collision curve (solid) in (t, a)-plane of Fig- 
ure 4, projected into the t) -space. The value of a is coded by the gray value. 
Values of other parameters: e = 0.1, = —0.1. The family of invariant curves is the 
transparent cone. The dashed curve inside the cone are the unstable equilibria of F. The 
fanned out intransparent surface are the switching lines. The points where the invariant 
curves touch the switching lines are marked as a solid black curve lying on the cone 
and the intransparent surface. 



Neimark-Sacker bifurcations and colliding smooth invariant curves can be efficiently computed 
directly because only the smooth two-dimensional map F and the expression for the boundary 
Do are involved. The Neimai^k-Sacker bifurcation in Figure 4 can be accurately and efficiently 
continued using standard numerical algorithms as described in detail in [20] and implemented in 
AUTO [13]. The smooth closed invariant curve is given imphcitly by the invariance equation 



>'('7(<P)) 



. T+?((p) 



y{(p), £ = [cosa,sina] •F_!*-'''V(<p) 



where y{(p) is given by 



and is the equilibrium of F 



yo 



y{(p)=yo + r{(p) 



cos<p 
sincp 



-Y^^'°yo, £ = [cos a, sin a] ■ Y!°yo. 



(22) 



(23) 



(24) 



The equations (22)-(24) have the periodic functions r{(p), ri{(p) and t{(p), the vector G and 
the scalar to as variables. The definition of y{(p) via (23) corresponds to a parametrization of the 
invariant curve by angle with respect to the equilibrium yo. The function r is the distance of the 
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point on the invariant curve from yo, rj is the circle map and t is the time elapsed from the last 
switch (the map F is defined only implicitly). This parametrization can be expected to work 
only in two-dimensional maps and for convex curves. Close to the point NSC in Figure 4 (in fact, 
close to the Neimark-Sacker curve) the invariant curves are ellipses and, thus, convex, making 
the parametrization (23) regular. We extend system (22)-(24) by the collision condition 

0= max t{(p). (25) 

(pe{o,27t] 

The right-hand-side of (25) is a smooth function of (p close to the point NSC in Figure 4 due to 
the convexity of the invariant curve. Thus, the system of equations (22)-(25) for the variables 
{r,ri,t,yo,to,z,a) defines a smooth family of colliding closed invariant curves which can be found 
by a Newton iteration with pseudo-arclength embedding [17]. Its projection into the (t, a)-plane 
is shown as a solid curve in Figure 4. The closed invariant curves of the colliding family projected 
into the space {x,x,r) are depicted (as a transparent gray cone) in Figure 7. Figure 7 gives 
evidence of the linear dependence of the radius r of the closed invariant curve on the parameter T. 
Thus, the solid collision curve must be quadratically tangent to the Neimark-Sacker curve at NSC 
in Figure 4 because the radius r has square root like asymptotics with respect to the distance of 
the parameter value to the Neimark-Sacker bifurcation. At the point ICB the invariant curve is 
close to a break-up such that the numerical approximation (32 complex Fourier modes) of r and 
T] gives an error estimate greater than 10^^. 

Remarks If the Neimark-Sacker bifurcation is subcritical then the line of collisions of closed 
invariant curves (solid) lies on the other side of the Neimark-Sacker bifurcation (hollow) in 
Figure 4 (but still tangential). The polygons remain stable filling the region between the solid and 
the dashed curve. 

The interaction between flip bifurcation and collision (dotted curves in Figure 3) has been 
analyzed in [ 1 8] . A periodic orbit L2 of period two branches off from L. In any generic two- 
parameter unfolding a curve of collisions of L2 emerges tangentially to the flip curve of L (looking 
very similar to Figure 4, replacing the Neimark-Sacker curve by a flip curve and the collision of 
the invariant curve by a collision of L2). The dynamics of the system between the collision of 
L2 and the collision of the period-one orbit L has, for example in [18], a stable period-two orbit 
flipping between the subdomains D and D+. 

7 Conclusions and open problems 

The local bifurcation theory of periodic orbits in hybrid dynamical systems with delayed switching 
can be reduced to the study of low-dimensional maps because the local return map (Poincare 
map) of periodic orbits is finite-dimensional. The codimension one events for generic systems 
(corner collision and tangential grazing) were studied in [27]. This paper studies the case of 
simultaneous corner collision due to reflection symmetry, a case that has been encountered 
frequently in example studies [4, 5, 7, 15]. 

The presence of hysteresis in the switch simplifies the proof of the reduction theorem (Theo- 
rem 9), given here for the first time. However, the statement of the reduction theorem remains 
valid also without hysteresis. 
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The derivation of Theorem 9 can be generalized in a straightforward manner to other discrete 
symmetries and other than binary switches (with more notational overhead) as long as the 
symmetry can be reduced locally near the periodic orbit. 

The reduction extends the appUcability of theory and numerical methods that have been 
developed for smooth and piecewise smooth low-dimensional maps to systems with delayed 
switches. On the numerical side this includes direct continuation of periodic orbits and their 
bifurcations and discontinuity induced events (such as grazing and collision) and the continuation 
of smooth invariant curves. Robust and universal methods for continuation and detection of 
discontinuity induced bifurcations for periodic orbits have been developed by Piiroinen [24, 25]. 
Methods for the continuation of closed invariant curves have a longer history. Schilder [26] and 
Thakur [8] present recent implementations, the papers also give surveys of earlier work. There 
is, however, a large gap between recent developments of numerical methods for closed invariant 
curves and piecewise smooth systems and their actual availability in the form of software. Due to 
this gap the investigation of the oscillator in Section 6 could not rely on generally available tools. 

The analysis of the oscillator shows that dynamical phenomena of hybrid systems with delayed 
switches can be systematically discovered with the help of numerical continuation and the 
reduction theorem. Two major open problems can be identified from the results of this analysis 
and the comparison to other results. 

First, there is a gap between the relatively simple local bifurcation theory of periodic orbits as 
presented here and in [27] and the abundance and variety of complex phenomena observed in 
detailed example studies such as [4, 7, 15]. In particular, many of the complex phenomena (say, 
chaotic attractors that have periodic orbits with various switching patterns embedded) cannot 
be reached by systematic continuation of periodic orbits and branching off at continuous local 
bifurcations. One of the reasons behind this gap is the presence of discontinuous events such as 
the grazing and the corner collision with ^ < presented in Section 4. 1 . These events are beyond 
the scope of local bifurcation theory because some trajectories leave the local neighborhood 
of the periodic orbit. A way to close this gap may be the study of discontinuous events under 
the assumption of the existence of a smooth 'global' return map for trajectories leaving the 
neighborhood similar to the treatment of global bifurcations of smooth dynamical systems. 

The second open question is the robustness of discontinuity induced bifurcations with respect 
to singular perturbations. For example, do the closed invariant polygons (as shown in Figure 5) 
persist when we change the oscillator (1) to 

x+^y + (0^x = u, I5y = x — y (26) 

(a modification to non- viscous damping for visco-elastic materials [1]) with a small j8? In 
particular, are there still finitely many smooth arcs? Note that the argument of [31] and of 
Section 6 that the arcs are images under F of the one-dimensional image rgF+ cannot be applied 
anymore because dim TgF+ = 2 for (26). Most statements of the bifurcation theory for piecewise 
smooth systems cannot be easily generalized to higher dimensions due to the lack of center 
manifolds. For example, the limit of small delay T, which is governed by standard singular 
perturbation theory for smooth dynamical systems [6], is highly non-regular for (1) in the case of 
£ = 0, see [27]. 
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A Basic properties of the forward evolution 

Let (^0 £ C([— 0,O];M") be the initial history segment of the continuous variable y (where > T) 
and Mo G { — 1, 1} be the initial state of the discrete variable u in system (4), (5). How does this 
initial state evolve to time T (defining the evolution E^ {t,{),uo)l 

First, we define £'^(,^o,mo) for T E (0,t] using the variation of constants formulation of (4). 
We define the following subsets of the closed interval [— T,0]: 



The set Rq is open relative to [— T,0]. Thus, if non-empty it is a union of countably many disjoint 
open (relative to [— T,0]) intervals. We arrange this sequence of countably many disjoint intervals 
into two subsequences of intervals: if {j = I, . . .), ly (j = I, . . .) and, possibly, one extra interval 



R+ = {se [ 
R- = {sG[ 
Ro ={s£ [ 



r,0]:h{^o{s))<-e} 
z,0]:h{^o{s))>e} 



)}■ 



ij' = intervals of Rq that have a lower boundary £ 7?+, i.e., h{^o{sJ)) = —e. 



Ij = intervals of Rq that have a lower boundary Sj G R , i.e., h{^o{Sj )) = e. 




( [— T,?o) if — T G Rq (then to is the upper boundary of the left-most interval of 
k = I Ro) 
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We define the following function ;^ : [0, t] M: 





' 1 


if 5- 






-1 


if 5- 


T£R-U[jjIJ 




, "0 


ifs- 





(27) 



Thus, X is measurable on [0,t] and either 1 or —1 everywhere. The variation-of-constants 
formulation of (4) is 

y{t) = U0) + -1) ~f (28) 

This is a fixed point problem on the space C( [0, t] ; M") of continuous functions on the interval 
[0, t] that has a globally unique solution y{-) G C([0, t] ;M") due to the Lipschitz continuity of 
/(•,±1). Then, for r G (0,t], £^((^o,mo) = {^{T){-),u{T)) is defined by 



yiT + s) ifse[-T,0], 

yo{T+s) ifse[-T,-T], (29) 



We observe that the initial value of the discrete variable uq only affects the result E^{yo,UQ) if 
/j(i§o(— t)) G (—£,£). Otherwise, (29) simply sets the discrete variable to its consistent value. 
For r > T we define E^{^o,uo) as a concatenation of smaller time steps, for example, if T G 
[{k—l)T,kT] thenf"^ =E'^ 1'^ o . . .oE'^ 1^ . This definition is independent of the particular partition 

of [o,r]. 

The definition of the evolution E using the variation-of-constants formulation (28) allows one 
to initialize E from arbitrary continuous history segments i^o and discrete states uq even if t,Q 
crosses the switching manifolds {h{x) = ±1} infinitely often or if mq is 'inconsistent'. Due to (28) 
the continuous component Ec of ^^(jOjWo) depends continuously on T. If T > T then E[{t,(),UQ) 
is Lipschitz continuous with respect to t. Its Lipschitz constant is Lmax • ll^cC'sO! "o) lre[ro,f£] II where 
^max is the Lipschitz constant of the right-hand- side /. For positive times the trajectory of the 
discrete component E^ is continuous from the right, that is limi.^,£'^(i^05 "o) = E'j{^o,uo) for all 
t>0. 



Proof of Lemma 1 (Point 1) Let 5 > be such that \s — t\ < 8 implies \h{E,Q{s)) — 
h{t,o{t))\ < le for all s,t G [-T,0]. This 5 exists because the initial value t,Q G C([-T,0];M") is 
uniformly continuous and h is Lipschitz continuous. 

We have to check how many sign changes the function defined in (27), can have in the 
interval [0, t]. Let us denote the upper boundary of each interval by f^. If /q is non-empty x 
can change its sign only in [fo, t] (one change is, possibly, in t^). If /q is empty we will use the 
notation ?o = in the following argument. 

After ?() the function x changes from —1 to 1 only at times t = tj +T when tj G R+, that is, 
h{t,f){s7j)) = £ (by definition of IJ) and h{^o{tj)) = —e. By definition of 5 this implies that 
tJ >sj + 5. Hence, X = ~^ for at least time 5 before it can switch to 1. The same argument 
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applies for intervals Ij and x switching from 1 to — 1 . Consequently, between two subsequent 
switchings of ;t; a time of at least 5 must elapse. This limits the number of switchings to a finite 
number on a bounded interval. 

Point 2: After time T the constant d in the above argument is bounded from below by 



Wmax-Lmaxmax||3;|| 

where y is the solution of the fixed point problem for the variation-of-constants formulation 
(28), //max is the Lipschitz constant of switching function h, and L^ax is the Lipschitz of the 
right-hand-side /. 

Proof of Lemma 4 (continuity) It is sufficient to prove the continuity of EJ{^,u) with 
respect to ^ in = y^{t + •) G C([— 0,O];M") and u = u^{t) for times T <r because (j'*,^*) is 
periodic. 

Due to the Lipschitz continuity of /(•, ±1) in the variation-of-constants formulation (28) it is 
sufficient to prove that, for any given A > 0, we can find a neighborhood such that 

r \E-^{^,u4t))-u4t + s)\ds <2A (31) 

JO 

for all 1^ E ?7 (i^*). That is, we have to show that, starting from ^, we follow the same flow as 
all the time in [0, T] except in a union of intervals of overall length A. Let t^ k(k= 1 , . . . , /i) be 
the crossing times of in [— T,0]. Let us denote by the value of the discrete variable at time 
t and at the switching times in [t,t + z], that is, uq = u^{t), Uj^ = u^{t + T + f* (^ = 1) • • • iM)- 
Since is continuous from the right u^{t + x + s) =Uk also for s slightly larger than — T and 
{k=\,...,n). 

Due to the weak transversality of (y* , ) for all sufficiently small 5i > there exists a ^2 > 
such that 

1. — ^2 > Uk-h{y^{t + t^ j^ + + e for all ^ = 1, . . . and 

2. — ^2 > uic-ih{y^,{t + s)) — e for all se [t* k-i,t* k — Si] (for ^ = 2, . . . ,/x) and s G [— T,f* i — 

Both statements follow from the strict monotonicity of |/z(3'*(-)| at crossing times of (y*,!/*) and 
the identities Myt/j(j*(f + + £ = and Uk-ih{y^{t + t^^k)) — ^ = 0. The points 1 and 2 imply 
that for all sufficiently small 5i > there exists an open neighborhood U {^^) such that £//(<§*) 
satisfy 

1. 0>Uk-h{^{u_k + Si)) + efora.\\k= l,...,;U,and 

2. > Uk-ih{^{s)) -E for all s £ [t*,k-\,t*,k- 5i] (for k = l,...,}X) and s £ [-T,f*,i -5i]. 

Consequently, for initial values ((§ , u) with ^ £U{^^) and u = u.^:{t) the discrete variable E^^{^,u) 
is equal to uq in [0,T + ?*j — di] and equal to Uk in [T+t^^k + Si,T + t^^k+i — 5i] for = 1, . . . ,jU — 1 
and, if t^^^ + 5i < 0, equal to in [r + f*,^ + 5i,t]. Thus, Ej{^,u) and u■^:{t + s) can be 
different only in + 1 intervals of length less than 25\. Choosing 5i = A/(/i + 1) we obtain the 
neighborhood U (<^*) guaranteeing (31) for all G ?7 
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Proof of Lemma 5 First we choose a 5 > sufficiently small such that |/j(y*(5))| is strictly 
monotone increasing in all intervals — 5,t^ + 5] where t^^t = !> • • • )M) the crossing 
times of the periodic orbit (j*,^*) in (— T,0). This 5 exists due to the weak transversality of 
{y^,u^). We define the constant 

co= min {\h{y{u,k - 5))\ - e\,\h{y{t.,^k + 5))\ - £\} , 

k=l,...,fl 

which is positive. Thus, \h{y^{t^^k — S) \ — e < —cq and |/j(3'*(?*,/t + 5))\ — £ > cq. We choose 
Ui such that for all ^eUi (C) 

max {\h{^{hM-d))-h{y4t^^k-d))\,\h{Ut*,k + S))-h{y{t,^k + d))\} <co/2. 

k=l,...,ll 

This guarantees that (•))!—£ changes its sign in — 8,t^ k + 5] for all G ?7 (i^*). Conse- 
quently, there exist unique times {k = \, . . . ,}JL) such that 

tk = mLn{t£[u^k-5,U^k + 8] ■ \h{^ {t))\ = e} , (k = l,...,n). 

If 5 is sufficiently small E'^{^,1) will change its value exactly at the times tk + T in [0, t] . This 
implies that = E^{^, 1) is given by the recursion (denoting to = — T and t^^i = 0) 

ryM*]^T(;^) iff G {tk,tk+i] and )t odd 
\Yl'"'^^^{tk) iff G {tk,tk+i] and k even. 

Hence, = ^^(i^, 1) depends only oni§(0) =3'(0) and fi, . . . ,f^. The discrete variable £^(1^ , 1) 
equals u^{t). Therefore, £''(1^, 1) = E^^'^ oE'^{^, 1) also depends only on ^{0),ti, . . . ,t^. 

Denote the continuous component Ec{£, , 1) of £'''(1^ 1 1) by The cross section £ is transver- 
sal to at 3'*(0). Thus, for every y € U (3'*(0)) there exists a locally unique traveling time t 
such that G £. The function tz'-y^tis well defined (and smooth) in ?7(y*(0)). If ?7i(<§*) 
is sufficiently small then ^''{0)^11 (3'*(0)) due to the continuity of Ec in = y.^,{-). The return 
time T{^) is given by T{^) = p + tz{^P{0)). As depends only on (i^(0),fi, . . .,t^) the same 
applies to T{^) and, hence, the return map P. 

Proof of Corollary 6 The condition 1 of Corollary 6 guarantees that the periodic orbit follows 
one of the flows in each of its crossing times f^;t G {0,p) (k = l,...,m), say, (Yk = Y±). 
Condition 2 guarantees that the flow F/t intersects the switching manifold {\h\ = £} transversally 
iny*{h,k) for all crossing times f* j., that is, h' {y^{t^^k))Yl{y*{t^.^k))\t=o / 0. 

The coordinates yo G L and (f 1 , . . . , ) of an initial condition in , 1 ) G o5^i must be in a small 
neighborhood of {y*{0),t^^m-^i+i — P, - ■ ■ ,f*.m — p)- Therefore, the headpoint y{t) of E[{t, , 1) 
follows Yk at times near t^ ^ for k = I, . . . ,m (because y^ follows yjt). Thus, the transversal 
intersection of Y^ with the switching manifold near f* 1, . . . ,„ implies that the crossings of 
y{t) near , . . . ,f*.m depend smoothly on the coordinates (yo,^!, • • • i^/x)- Consequently, also 
all times t G {0,p) when E'^{t,,\) changes its value depend smoothly on {yQ,t\,...,t^). As 
the cross-section Z is transversal to Y+ and y{t) follows Y+ for t near p the return time also 
depends smoothly on {ya,t\, ... ,t^) and, hence, the whole return map P depends smoothly on the 
coordinates jo £ ^ and (f 1 , . . . , ) . 



^(0 = 
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B Proof of Lemma 8 and Theorem 9 



Proof of Lemma 8 The strict transversality condition (13) implies that the periodic orbit also 
satisfies the weak transversality condition (as given in Definition 3). In the proof of continuity 
(Lemma 4) we established that the discrete components E[^{t,,u) of trajectories starting from 
initial conditions {t,,u) near ()'*(•) )"*(0)) change their values always close to times where 
changes its value. Also the direction of change must be the same because there is a minimum 
distance between subsequent changes (given in (30)). The colliding symmetric periodic orbit 
changes its value exactly twice per period, at f = and at f = T. Thus, the image {t, , u) of an 
initial value (i^o, "o) must have a continuous component t, switching exactly once near s = —A 
from Y to F+. (The Poincare section £ was taken at 3'*(A) and switches from —1 to 1 at 
f = 0.) 

Consequently, must have the form (14) for some y € ?7 (3'*(0)) and a time 6{y), which is the 
traveling time from y to £ following F+. The only open question is if i§ (0) G £ and the time t\ ) 
are uniquely determined by y. 

Let y ^U{y^ (0) ) be given. The time d {y) is implicitly given by Y^y G Z, which means 

fI[Yty-y.m=0. (32) 

The linearization of (32) my = y^ (0) with respect to 6 is f2 fi, which is non-zero due to the strict 
transversality (13). Thus, is a locally unique and well defined function of y. Consequently, 
t, (0) = Y^y is also well defined and smoothly dependent on y. The time ?i = ?i (i^ ) is given by 
h =t{y) — Q{y) where t is implicitly defined by 

£ = h{Yly) ifh{y)<e, (33) 
e = h{Yly) ifh{y)>e. (34) 

The linearization with respect tot my = y*(0) is in case (33) and h'^ f\ in case (34). Both 
linearizations are non-zero due to the strict transversality condition (13). 

Proof of Theorem 9 Let y (3'*(0)) be given. Lemma 8 gives a unique element of the 
image of P corresponding to y that has the form (14). The image of y under the map m is the 
location y2 at the next time 52 when the discrete component E^^{t, , 1) changes its value from —1 
to + 1 . The full reflection symmetry of the periodic orbit L implies that m is the second iterate of 
a map F •.U {y^{Q)) ^ U (3'*(0)) which is defined as —y\ where y\ is the location at the next time 
si when the discrete component £'^' (i^ , 1) changes its value from +1 to —1. 

The regular implicit condition (32) defines d{y) and §(0). The regular implicit condition (16) 
(same as (33), (34)) defines the time ti{y) = t{y) — d{y) locally uniquely. Thus, the 1) 
changes its value to —1 at time ^'i = x + t\ {y). The point y\ is the headpoint of the continuous 
component Ep{t,,\). It has the form 

3;i = F|i§(0) = F|'+% = F^+'Wy 

because ?i = t{y) — Q{y) where t{y) is given by the regular implicit condition (16). Thus, 
F{y) = —y\ has the form claimed in Theorem 9. 
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